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RESONANT EQUILIBRIUM CONFIGURATIONS IN 
QUASI-PERIODIC MEDIA: PERTURBATIVE EXPANSIONS 

RAFAEL DE LA LEAVE, XIEENG SU, AND LEI ZHANG 


Abstract. We consider 1-D quasi-periodic Erenkel-Kontorova models. 

We study the existence of equilibria whose frequency (i.e. the inverse 
of the density of deposited material) is resonant with the frequencies of 
the substratum. 

We study perturbation theory for small potential. We show that there 
are perturbative expansions to all orders for the quasi-periodic equilibria 
with resonant frequencies. Under very general conditions, we show that 
there are at least two such perturbative expansions for equilibria for small 
values of the parameter. 

We also develop a dynamical interpretation of the equilibria in these 
quasi-periodic media. We show that equilibria are orbits of a dynami¬ 
cal system which has very unusual properties. We obtain results on the 
Lyapunov exponents of the dynamical systems, i.e. the phonon gap of 
the resonant quasi-periodic equilibria. We show that the equilibria can 
be pinned even if the gap is zero. 

Quasi-periodic Frenkel-Kontorova models, resonant frequencies, equi¬ 
libria, quasicrystals, Lindstedt series, counterterms 

[2010] 70K43, 37J50, 37J40, 52C23 

1. Introduction 

The goal of this paper is to formulate the theory of resonant equilibria 
in quasi-periodic Frenkel-Kontorova models. We argue that these equilibria 
play an important role in the phenomenon of pinning (the equilibria that 
survive after an external force is applied). In periodic Frenkel-Kontorova 
models, the role of resonances in pinning is well established. 

We recall that in Frenkel-Kontorova models [iBK041ISel92|| , one considers 
configurations given by a sequence of real numbers (think of the position of 
a sequence of particles deposited on a 1-D material). The (formal) energy of 
the system is the sum of a term of interaction between nearest neighbors of 
the deposited material and a term modeling interaction with the media. In 
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the quasi-periodic Frenkel-Kontorova models studied here, the interacting 
potential will be a quasi-periodic function of the position reflecting that the 
medium is quasi-periodic. We will be interested in equilibria, i.e., config¬ 
urations such that the derivatives of the (formal) energy with respect to the 
position of each of the particles vanish. We note that even if the energy is a 
formal sum, the equilibrium equations are well defined. More details of the 
models will be discussed in Section 

In ||SdlL12b[ ISdlL12al , one can find a rigorous mathematical theory of 
quasi-periodic solutions whose frequency is not resonant (indeed Diophan- 
tine) with the frequencies of the substratum. The rigorous theory of IISdlL12bi 
|SdlL12al also leads to efficient algorithms that can compute these quasi- 
periodic solutions arbitrarily close to their breakdown. Implementations of 
these algorithms and investigation of the phenomena at breakdown appear 
in HBdlLlSH . The paper IISdlL12bl . studies models with nearest neighbor 
interaction while [|SdlL12all studies the case of long range interactions. 

Our motivation is to study the phenomena of “depinning”. When we add 
an external force (no matter how small) to the model, many quasi-periodic 
solutions disappear. Nevertheless, there are still other quasi-periodic so¬ 
lutions survive. This physically corresponds to the deposited material rear¬ 
ranges itself to withstand the force. It is well known in the periodic Frenkel- 
Kontorova models that the solutions which persist under forcing are reso¬ 
nant with the media. The papers IISdlL12bl ISdlL12all also show that, in the 
quasi-periodic case, the smooth non-resonant solutions do not exist when 
there is an external force. Hence we are interested in resonant solutions in 
the quasi-periodic models. 

In this paper, we show that there exist at least two formal power series 
in the amplitude of the coupling describing resonant solutions even in the 
presence of external forces. The delicate analytic question of convergence 
of these series will be pursued in HdlLSZlSH . We also develop a dynamical 
interpretation of the equilibria as orbits of a dynamical system. Using this 
dynamical interpretation, we show that the phonon gap of quasi-periodic 
equilibria is zero. 

We hope that this paper can lay the ground work for future explorations. 
We will use rigorous mathematical tools to study the convergence of the 
formal series developed here in HdlLSZlSl . We will also use numerical 
methods to explore in a non-rigorous but more quantative way some of the 
phenomena discussed here. 

Equilibria in quasi-periodic media with a resonant frequency have been 
investigated numerically in nvEFRJ99ilvEFJ01[|vEF02ll . These papers also 
studied the phonon gap and found it to vanish when there are smooth solu¬ 
tions (in agreement with the results here). 
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The variational and topological methods that have been proved useful in 
the periodic case do not extend to the quasi-periodic case in its full strength. 
Several interesting counterexamples are in ULSOSi IFed75H . Partial results 
are in [IGGP061lAPT0llGPT131[EOT3l . 

2. Models considered and formulation of the prob l em 

We consider models of deposition in a quasi-periodic one-dimensional 
medium. 

If An 6 R denotes the position of the n-th particle of the deposited mate¬ 
rial, the state of the system is specified by a configuration, i.e. a sequence 
{^Jnez- We associate the following formal energy to a configuration of the 
system 

(1) J^(x) = ^ ^(a„+i - Xn - af - V(Xna) - AXn 

neZ 

where V : T'^ ^ R is an analytic function, or 6 R'^ is an irrational vector and 
a, d are some real numbers. 

A natural example often used is 


(2) ^(a) 


- ^(a„+i - x„ - a)^ + a cos x„ + B cos( V2a„ + (f>) - 


AXn 


which corresponds to cr = (1, V2) and V(0i, 62 ) = A cos 61 + B cos(02 + 0)- 

The term (a„+i - x„ - af represents the interaction among neighboring 
deposited atoms. The term V{x„a) represents the interaction with the sub¬ 
stratum. The interaction at position a 6 R is the quasi-periodic function 
y(AQ;). This models that the substratum is quasi-periodic. The term Axn has 
the interpretation of a constant field applied to the model. In the case of 
deposited materials, we can imagine that the sample is tilted and A is the 
component of the gravity. The physical meaning of the equilibria corre¬ 
sponding to d ^ 0 is that, when we apply a small external field, the config¬ 
urations rearrange themselves so that they can respond to the force and do 
not slide. Of course, when the force is strong enough, no rearrangement is 
possible and all the configurations slide. This is the microscopic origin of 
static friction. 

In the case of periodic media, the pinned solutions are known to corre¬ 
spond to resonant frequencies. Since the resonant tori do not survive, it is 
reasonable to consider resonant frequencies. 

The existence of external forces d is a very important novelty with re¬ 
spect to the previous papers [ISdlL12bilSdlL12aL It was shown in [ISdlLI2bl 
ISdlLI2al that if there is non-resonant quasi-periodic solution, then d = 0. 
In our case, we will show how to construct quasi-periodic equilibria with 
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nontrivial A and will show how to compute perturbatively the range of sueh 
A for whieh solutions with a preseribed resonant frequeney exist. 

Without any loss of generality, we ean assume that 

(3) k-aiN Vyt6Z'^-{0}. 

If there existed a resonanee k ■ a = 0, we eould just use less frequeneies to 
express the quasi-periodie funetion. 

2.1. Equilibrium equations. A eonfiguration is in equilibrium if the forees 
aeting on all the partieles vanish. Equivalently, the derivatives of the energy 
with respeet to the position of the partieles vanish. That is, 

—— (x) = 0 V n 6 Z. 

OXn 

In the model O, the equilibrium equations are 

(4) x„+i + x„-i - 2x„ + daV(x„a) + 4 = 0 VneZ 

where da = a - V and V is the usual gradient. 

Note that even if the energy O is just a formal sum, the equilibrium 
equations © are well defined equations. 

It is very tempting to eonsider © as a dynamieal system, so that we 
obtain Xn+i as a funetion of and x„_i. This system has very unusual 
properties. This will be pursued in Seetion[5l 

2.2. Quasi-periodic configurations, hull functions. In this paper, we will 
be interested in quasi-periodie solutions of frequeney m 6 R. 

These are eonfigurations of the form 

(5) x„ = nco + h(n(oa), 
where h : T‘^ ^ R. 

A eonfiguration given by a hull funetion © satisfies the equilibrium 
equation © if and only if the hull funetion h satisfies 

(6) h{n(oa + coa) + h{na)a - coa) - Ihincoa) -l- daV{ncoa + ah(na)a)) -l- d = 0. 

The equation ® was eonsidered in [ISdlL12bllSdlL12all when oja is Dio- 
phantine (in partieular, ncoa is dense in the torus T^^). 

In our ease, ncoa will not be dense on the J-dimensional torus (see See- 
tion 12.31) and the equilibrium equations we will derive are different from 
those in [|SdlL12bl |SdlL12a| . 
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2.3. Resonances. The goal of this paper is to study situations when there 
are k - {0} and m 6 Z sueh that 

(7) k ■ coa - m = 0. 

When dV]) holds we say that (k, m) is a diserete resonanee for coa and we 
refer to the pair {k, m) as a resonanee. 

Remark 1. Note that these discrete resonances (|7]) are different from the 
resonances of the media we excluded before (k ■ a Q, 'i k e - {Q}). 

Remark 2. If 

k-ai^O MkeZ^^X {0}, 

given any kg eZ‘‘ \ {0}, m eZwe have that co = -ml{ko ■ a) is a resonant 
frequency. Since k^-a can be arbitrarily large, we see that the set of resonant 
frequencies is dense on the real line. Of course, once we fix a, the set of 
resonant co is a countable set. 

2.3.1. Multiplicity of a resonance. Clearly, if (k, m), (k, in) are diserete res- 
onanees so is {k + k,m + m). 

In mathematieal language, 

‘^oja = {(^5 m) eZ‘^ xZ : k ■ coa - m = o| 

is a Z-module ealled the resonanee module for co. 

We denote by l{co) = dim{.y^^,Ja) the dimension of the resonanee module 
and we eall it the multiplieity of the resonanee. The meaning of l{co) is 
the number of independent resonanees. We ean find (/ci, mi),..., (fc/, mi) in 
sueh a way that all resonanees ean be expressed as eombinations of the basie 
resonanees (and also no other set of basie resonanees with smaller number 
of elements will allow to express all the resonanees). 

2.3.2. Only resonances of multiplicity 1 appear in the models dU). In Hamil¬ 
tonian meehanies for systems with d degrees of freedom, one ean find res¬ 
onanees of all multiplieities up to d. As we will see later, in Seetion [51 
one ean give a dynamieal interpretation of the equilibrium equations as a 
dynamieal system in J -l- 1 dimensions. Nevertheless, in our models only 
/ = 1 appears independently of the number of degrees of freedom. This 
highlights that the problem here is different from the Hamiltonian problem. 

Proposition 1. If coa is resonant, i.e. + {0}, then l(co) = 1. 

Proof. Note that 

ki • coa - mi = k 2 ■ coa - m 2 = 0 
implies (beeause mi 0, m 2 0 beeause of dS])) 

mi m 2 

ki ■ a k 2 ■ a 


CO = 
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and therefore 


a ■ (kim2 - k2tni) = 0 

and, beeause a is non-resonant ([3]) we have 


kim2 = ^ 2 ^ 1 - 


Therefore, the two resonant veetors are related. □ 

2.3.3. The intrinsic frequencies. When coa is resonant, we ean find a matrix 
5 6 S L{d, Z), 6 , L 6 Z^ in sueh a way that 

( 8 ) Bcoa = {Q, 0) + L with kl-kiZforke Z^“' - {0}. 

We will refer to f2’s as the intrinsie frequeneies. They are essentially 
unique, i.e., unique up to ehanges of basis in R^“^ given by a matrix in 
SL{d-1,Z). 

In this ease, the set {na)a]„e'z has a elosure whieh is a J - 1 dimensional 
torus. This torus is invariant under the translation Ta)a- If we stay in this d- 1 
dimensional torus, T^a ean be deseribed as Tq. The torus T"' is foliated by 
these indexed by another parameter 77 6 We will write a point in 
as {if/, 77 ) where if/ is the eoordinate eorresponding to the position in 
The eoordinate 77 seleets the d - \ torus we are eonsidering. 

2.4. Quasi-periodic equilibria with resonant frequencies. The natural 
notion of the hull funetions in the resonant ease would be to assume that the 
equilibrium solutions have the form 

(9) Xn = noj + v{nkl) 

with V : R. 

Note that the physieal meaning of co is still the mean spaeing of the solu¬ 
tions (i.e., an inverse density). The term v{nQ.) represents fluetuations that 
ean be parameterized in terms of the intrinsie frequeney Q.. Of eourse, we 
eould represent them in terms of the original frequeneies, but it is more 
natural to ehange variables so that they beeome a part of the equation. 

We will refer to the 77 variable as the transversal phase. The resonant 
solutions eonsidered here, eover densely a torus of eodimension one. The 
one-dimensional variable 77 measures the position of these eodimension-one 
tori on the eonfiguration spaee eorresponding to the internal phases of 

y. 

If we substitute the parameterization ® into the equilibrium equation, 
we obtain that the equilibrium equation dH) is equivalent to: 


(10) v(nQ. -I- Q) -I- v{nkl — Q) — 2v(nQ) -I- daVincoa -l- av{nkl)) -l- d = 0. 
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If we furthermore introduee the notation daV(6) = W(B6) and Ba = 
yS, and observe that the nO. is dense on T‘^“\ we see that for eontinuous 
funetions v, (fTOl) is equivalent to: 

(11) v{i/r + Q.) + v{i^ - Q.) - 2v((A) + 0) + ySv((A)) + d = 0. 

Note that we ean also eonsider solutions of the form 

(12) x„ = no) + v(nQ. + ^i) + ^2 

for any fixed and ^ 2 - This will give us freedom to add the transversal 
phase 77 as an additional parameter in (fTTI) . So the equilibrium equation we 
will eonsider is 

(13) v(i/^ + Q) + v(i/r - Q) - 2v(i/f) + W((i/f, rf) + ySv(i/^)) + d = 0. 

By simple ealeulations, it ean be shown that (fT3]) is equivalent to dH) when 
the hull funetion is of the form (fT^ and Ba ^2 = (^ 1 , 17 ) ■ 

Remark 3. Because P has components both in the if/ and the rj directions, 
the equation (fT3l) cannot be considered as a parameterized version of the 
equations considered in nSdlL12bl . As we will see, the symmetries of the 
equation involve transformations that mix the dependence in if/ and in q. 

2.5. The symmetries of the invariance equation (fT3]). The equation (fT3]) 
possesses remarkable symmetries that make the solutions not unique. These 
symmetries lead to Ward identities. In eontrast with the ease of non-resonant 
solutions, the group of symmetries is infinite dimensional. In [ldlL08[[S^dlL12bl 
ISdlL12al these symmetries are used to develop a KAM method. 

The main observation is that if (v, d) is a solution of (fT3]) . then, for every 
dq) : ^ R, the pair (v, d) is also a solution of (fT3]) where we denote 

p = and v, d are defined by: 


v{iff, q) = v(((A, q) + dqffS) + L{q), 
Xq) = d(77 -P L{q)ff}. 


(14) 


Notiee that the symmetry (fT4l) involves ehanging not only the argument 
if/ but also the argument q. Note the spaee of symmetries of the equation is 
not just a finite dimensional spaee but rather an infinite dimensional spaee 
of funetions. 

2.6. A normalization of the solutions of the invariance equation (fT3l) . 

For later applieations, it will be useful to have loeal uniqueness of the so¬ 
lutions (e.g. to diseuss smooth dependenee on parameters, perturbative ex¬ 
pansions on parameters). Henee we impose the normalization 


( 15 ) 
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Since the symmetry (fT4l) involves ehanges of arguments, giving a v,,, 
finding the L{r]) that aeeomplishes the normalization involves solving the 
implieit equation 


(16) 


/(?7 +/?;,t( 77 )) + L{rf) = 0 


where lij]) = v((A, rj) dif/. 

If I and its derivative are small, one ean solve eqrefimplieit using implieit 
funetion theorem. 

2.7. Diophantine condition. In eontrast with KAM theory, we will not 
need very delieate estimates on the solutions and henee, we ean deal with 
very general Diophantine eonditions. We will assume that Q. satisfies 


(17) 


lim — sup \n\k ■ Q. - m\ = 0 . 

N-^OO N \k\<N,meZ 


Note that the eondition (fT71) is mueh weaker than the usual Diophantine eon¬ 
ditions and even than the Bjruno-Riissmann eonditions. The eonditionlfTTl) 
is the natural eondition in the study of existenee of series to all orders. 

The following proposition shows that the sets of frequeneies we are eon- 
sidering are abundant. 

Fix a veetor k 6 \ {0}, m 6 Z \ {0} and assume without loss of generality 

that there is no eommon divisor in the eomponents of k. For any a 6 
satisfying a ■ k 0, we ean find a unique co sueh that a ■ koj - m = 0. 
Fix a Bk & S L(d, Z) sueh that the last row of B equals k (It’s possible to 
find sueh B sinee the eomponents of k has no eommon divisor). Then, let 
Bkaoj = (D, 0) -I- m. Henee, for any k, m and a satisfying a ■ k 0, we ean 
define Q as a funetion of a. Denote Q = Fk,m(o:). 

Proposition 2. The set of a for which Fk^,„(a) satisfies (fT71) for all k,m is 
of full measure in R^. In other words, for a full measure set of medium 
frequencies, we can find a countable many resonant frequencies that lead 
to intrinsic frequencies satisfying (fTTI) . 

Proof Sinee eountable interseetions of sets of full measure are of full mea¬ 
sure, to prove Proposition [2] it suffiees to show that for a fixed k, m as above, 
the set {or e R'^la ■ k 0 and Fk,m{(F) satisfies (fT7]) l is of full measure. 

Beeause the set of D’s whieh satisfy (fTTI) is of full measure on R"^"^ and 
the linear map Bk is differentiable and surjeetive, the preimage of the set 
of D’s that satisfy (fTTI) under Bk is also of full measure in the hyperplane 
F = {y|y-k-m = 0}. Denote the set of preimages as F'. Then any nonzero 
sealing of an element of F' will give an a we want, whieh also form a full 
measure set in R'^. □ 
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3. Function spaces and linear estimates 

The main tool that we will use to eonstruet perturbation theories is the 
solution of eohomology equations. 

We denote 

Dp = {0 6 C‘l7f I |lm(0,)l < p] 

and denote the Fourier expansion of a periodie mapping rf) on Dp by 

keZd 

where • is the Euelidean sealar produet in and Vk are the Fourier eoeffi- 
eients. 

We denote by the Banaeh spaee of analytie funetions on Dp whieh 
are real for real argument and extend continuously to Dp. We make s^/p a 
Banaeh space by endowing it with the supremum norm: 

llvllp = sup_|v(iA,77)|. 

These Banaeh spaees of analytie funetions are the same spaees as in 
[|Mos67ll . 

We will eonsider equations of the form 

( 18 ) viil/ + n,T])- 77) = 0((/r, ?7), 

where if/ 6 

To simplify our notations, we will denote v(i/^ + £l) and v{if/ - Q.) as v+ 
and V-, respeetively. Similar notations will be used for other funetions. We 
also use T to represent the translation operators, i.e., Tav{i]/) = v(i/^ + Q.). 

Lemma 1. Let (p g be such that 

(19) r = 0, 

for all rj. 

Assume that Q. satisfies the assumption (fT71) . 

Then, for a fixed ij, there exists a unique solution Vrj of (fT^ which satisfies 

(20) r v{fi, q)dil/ = 0. 

The solution v 6 s^p> for any p' < p and we have 

llvllp/ < C(p,p')\\fp\\p. 

Furthermore, any distribution solution of (fTSl) differs from the solution 
claimed before by a constant. 

Iff is such that it takes real values for real arguments, so does v. 
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If we consider now the dependence in tj, we have that v e and 

llvllp/ < C(p,p')||0||p. 

Proof. We note that, as it is well known that obtaining v solving (fTSl) for 
given (f) is very explieit in terms of Fourier eoeffieients. If 


then, V is given by 


v(«A, 77 ) = ^ 

ki^^ ki^0,m 


irv 


2m{k-il/+mr]) 


Using Cauchy estimates for the Fourier coefficients \fk,m\ ^ exp(-2:?rp(|k| + 
Nl))ll0llp and that - l|-i < Cdist(k • Q, Z)“^ and the assumption (fT^ . 
we obtain that 

llvllp. < C ^ exp(-27rp(|k| + |m|))||0||pdist(k • Q, 

ki=Q,m 

< C\\(p\\p ^ exp(-2;rp(|k| + |m|))dist(k • Q, Z)“^ exp(2;rp'(|k| + |m|)) 

ki=Q,m 

< c(p,p'm\p. 


□ 


4. Lindstedt series eor quasi-periodic solutions with resonant 

FREQUENCIES 

The goal of this section is to study (fTSl) perturbatively when the non¬ 
linear term is small. Hence, we will write (fT3l) with a small parameter s 

( 21 ) 

vif + Q.,ri) + vif - Q, 77 ) - 2 v(i/ 7 , 77) -F eWHf, 77) -F 77)) -F A{p) = 0. 

We will find v{f,p),A{p) solving (l 2 TI) and (fTSl) in the sense of formal 
power series in e. In this paper, we will not consider the problem of whether 
these series converge or represent a function. This will be studied in more 
details in HdlLSZlSH . 

Since one may want to find solutions correspond to d = 0 (or d = d* with 
|d*| small), it is important for us to keep track of |^( 77 , e) in order to solve 
d( 77 , 6 ) = 0 by implicit function theorem. 
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Following the standard perturbative proeedure we will write 


oo 


( 22 ) 



oo 


oo 


A=y 6 "^”. 


n=0 


Here v" and A" are eoeffieients of e", not powers of v or A. Substitute (l22l) 
in (I 2 TI) and equate powers of e. 

Of eourse, earrying out this proeedure for n < N will require that Q. 
satisfies some Diophantine properties as well as some differentiability as¬ 
sumptions. 

Equating the eoeffieients of 6 ** in (|2TI) we obtain 


(23) v°(i/r + Q.,r]) + -Q.,ri) - 2v°(^, 77 ) -l- A^ii]) = 0. 


Henee, if Q satisfies the eondition (fTTI) we see that is eonstant, /l° = 0 
and imposing the normalization (fT5l) we obtain v° = 0 . 

Matehing eoeffieients of in both sides of (l2TI) we obtain 

(24) v^(i/r + Q,jj) + v\t/f - Q,7]) - rf) + VF(i/^, rf) + /l'(?7) = 0. 

We see that, using the theory in Seetion[3l to have analytie solving (l24l) . 
it is neeessary and suffieient to have 

(25) A\ri) = - f 

Then, v\ A^ ean be determined uniquely up to a eonstant from (l24l) . In 
faet, in Fourier series, the equation for v^,A^ is 


v^ 2 (eos( 2 ;rm) - \) =-Wk - 6o,kA\ 


(26) 


where So^k is the Kroneeker delta. In partieular, the eonstant in v* is deter¬ 
mined by the normalization (fT5l) . 

Proeeeding to higher order follows the same pattern. We see that mateh¬ 
ing the terms of order e” in (|2T]) we obtain 


(27) v”(i/r + Q,r]) + v”(i/r -n,7])- 2v"(i/r, 77) -p 77 ) -r A"ir]) = 0, 


where R” is a polynomial expression in ..., with eoeffieients whieh 
are derivatives with respeet to if/ of rf) + /3v{il/, 77 )). This polynomial 
ean be eomputed explieitly beeause it is given by 


( 28 ) 
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and these are well known formulae. We also note that, from the algorith- 
mie point of view there are effieient ways to eompute using methods of 
“automatie differentiation” HHarl ll[BCH^06L 

Sinee i?” ean be eomputed explieitly, (1271) ean be solved the same way as 

dal). 

We have therefore established 

Theorem 1. Assume that Q satisfies (fTTI) and that W : Dp ^ C is an 
analytic function. We can find formal power series solutions in e of the 
form (1771) solving the equation (I7TI) . For any 0 < S < p, each of the terms 
q) is analytic in Dp_g. If W takes real values for real values, then so 
do v”, and A is real. 


4.1. The auxiliary equation. Now, we turn to the problem of studying the 
equation 

(29) A{q,e) = A*. 

We expeet to obtain a solution q*{e) provided that (1791) satisfies some non- 
degeneraey eonditions. 

Having solution of (1791) to order 1 in e, amounts to 

A^(q) = 0. 


That is, we need to find q sueh that 


(30) 


f 


W({fi,q))dfi = 0. 


Theorem 2. The equation (l30l) has always two solutions. 


Proof. Sinee 

r W((fi,q))dfi = f (daV)(B~\fi, q))dfi, 

Jjd-I 

if we integrate again with respeet to q we obtain 

r r Wiiil/,q)dil/dq = f (daV)(B~f fi,q))dfidq = 0. 


(31) 


Henee the funetion of q given by q)dfi is a eontinuous periodie 

funetion of q with zero average. Therefore, it has at least two zeros. We 
also note that there are open sets of perturbations where there are 4,6,- • • 
zeros. □ 


Denote one of these solutions of (l30l) as q*. 

A suffieient eondition that ensures that we ean solve the equation (1791) to 
all orders is that 

d , I 

( 77 ? \7]=TJ*,€=0 


(32) 


^ 0 . 
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More explicitly, 

(33) r rj))dif/ dr] 0. 

Jjd orj 

Then, the implicit function theorem for power series [|Car95[ lDie711l gives 
us that we can indeed find 77 * (e). 

Similarly, we can solve the equation A{7]) = A* provided that |/l*| is suffi¬ 
ciently small. 

Therefore, we have established 

Theorem 3. Assume that Q e ^(v, r) as defined in (fT71) , that W is an 
analytic function, and that (l3^ holds, we can find formal power series rje 
in e so that is the solution of (l2TI) . 

Clearly, since the function A^lrf) are bounded, if A* - the physical force 
- is large enough, there is no solution. This has a clear physical meaning. 
If we increase the external force but keep it small, the system can react 
by changing the transversal phase. If the force increases beyond a thresh¬ 
old, the system cannot react by adapting the phase. Hence, the equilibrium 
breaks down. In this paper, we are not considering the dynamics of the 
model, only the equilibria (our models for the energy include only the po¬ 
tential energy of the configuration and not any kinetic energy). One can, 
however, expect that, if there was some dynamics, the equilibria considered 
here could slide. 

Of course, the sufficient condition (l3^ is far from being necessary and 
there are many other conditions that are enough. 

Proposition 3. Assume that Q satisfies (fT71) , that W{{tf/, rf) -l- rf)) is an 
analytic function, and that (13^ holds. 

Assume that rj* is such that for some m eN we have 

T(7/*) = 0, /=!,... 2m 

(34) 

A^'^^Ht]*) ^ 0 . 

Then, we can find formal power series in e so that v{f, pf) is the solu¬ 

tion of (I 2 TI) . 

The proof is again an application of the implicit function theorem for 
power series. 

5. A DYNAMICAL INTERPRETATION OF THE EQUILIBRIUM EQUATIONS OF 

Frenk ei, -Kontorova models 

In this section, we present a dynamical interpretation of the equilibrium 
equations dH) in Frenkel-Kontorova models. 
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Even if the dynamieal interpretation is possible for finite range inter- 
aetions, we see that adding another small interaetion of longer range is 
a singular perturbation (even the dimension of the phase spaee ehanges). 
Whereas, for the methods in this paper, adding a small term in the longer 
range is a regular perturbation of the same order. 

A straightforward way of transforming the equilibrium equation 

(35) x„+i + x„_i - 2x„ + edaViXnO) + /l = 0 'i n el. 


into a dynamieal system is setting 


(36) 


y„+i = {lyl -yl- €d„V{ayl) - X,yl). 


However, (l3^ is not very useful beeause we have to eonsider it as a map 
of and the term daV(ayl) does not make apparent that it is periodie in 

A more natural formulation is obtained by observing that the equation (jH) 
is equivalent to the system on x R 

Pn+l = Pn- edaViqn) - A 

(37) _ 


where qn 6 Pn 6 R. (Just multiply dH) by or and use the substitution 
Pn = Xn- x„-i, q„ = axn- Noto that (HI) is equivalent to 


(X„+1 - Xn) - {Xn - X„-i) + edaV{aXn) + d = 0 


henee, we obtain the first equation.) 

We will write the mapping (1371) as 

(38) {^Pfj-^.\, qfi^i) Fs^xi.Pn'^ 

Note that (1371) is typographieally very similar to the standard map IIChi791l 
or to analogues introdueed for volume preserving maps. Nevertheless, there 
are signifieant differenees (besides the different dimensions). 

A very erueial differenee between (l38l) and the generie volume preserving 
maps is that qn+i - qn is always a multiple of a (see (l371)l. So that the two 
dimensional leaves 


(39) Mq^ = {(p, qo + at)\p,te R} 

are preserved. Note that eaeh of the leaves Al^,, is dense in the d+\ dimen¬ 
sional phase spaee. 

The mapping (1371) elearly preserves the volume form dp /\ dqi A ... Aq^ 
sinee it is the eomposition of 

Pn+l = Pn- edaViqn) “ d 
QnA-\ ~ 


( 40 ) 
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and 


(41) 


Pn+l ~ Pn 

Cjn+l ~ ^Pn+\- 


We recall that, in our context, a volume preserving map is exact when 
F*{pdq\ A dq 2 A ... A dqd) = pdqi A dq 2 A ... A dqd + dP where P is d - I 
form. 

Indeed, (l38l) is an exact volume preserving map if and only if d = 0, since 
it is easy to observe that, when 4 = 0, both (l40l) and (1441) are exact. 

When e = 0, d = 0, the map (lJ7l) is integrable. The codimension-one tori 
given hy p = cte. are invariant and the motion in them is a rotation. 

The volume preserving KAM theory leads us to expect that for e 
1, d = 0, the tori in which the frequency of the motion is Diophantine sur¬ 
vive. We also expect that the tori with resonance, breaks down into lower 
dimensional tori. The lower dimensional tori can survive for |d| 1 (de¬ 

pending on e). 

In this paper we have quantitative (but formal and non-rigorous) predic¬ 
tion of these phenomenon based on perturbative expansions. We hope that 
some of them may be either verified by rigorous results or explored numer¬ 
ically. 

5.1. On the global geometry of the constraints given by (1^ . Integrable 
systems with constraints have been studied extensively in geometric me¬ 
chanics. Nevertheless, the systems we consider here have some unusual 
properties that we would like to highlight. 

It is customary to classify the constraints in holonomic when the distri¬ 
butions are integrable (in the sense that they foliate the phase space with 
a smooth quotient) and non-holonomic when the distributions are not in¬ 
tegrable and they violate the hypothesis of Frobenius Theorem IISou97[ 
IAud08[[HoTTT]l . 

The constraints (1^ escape this dichotomy. They are locally integrable 
(they do satisfy the hypothesis of Frobenius Theorem and are locally given 
by invariant manifolds that give rise to a foliation) but nevertheless, the 
manifolds are dense, so that they do not give a nice quotient manifold. 

Hence, even if we have holonomic constraints locally (and the infinites¬ 
imal results about holonomic systems are applicable), some global aspects 
such as symplectic reduction [|Mev73[ IMW74llMW01ll cannot be applied to 
(ETI). 


5.2. Lyapunov exponents and phonon localization. In this section we 
study the so called phonon gap around the equilibria of (fTOl) given by a hull 
function. 
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Let US Start by recalling some standard definitions. The main idea is that 
sound waves are defined by the propagation of infinitesimal disturbances 
around an equilibrium equation. 

If we linearize around an equilibrium solution x = {Xn}nez, we obtain the 
dynamics of the infinitesimal perturbations is given by 

(42) + (dJ^V(ax„)^„ = {LSn. 

It is clear that the propagation properties of sound waves will be affected 
by the spectral properties of the operator JLx- 

Note that the operator X is a one-dimensional Schrodinger operator with 
a position dependent potential. The dependence will be given by the dy¬ 
namics of the Xn- In particular, for the solutions given by a hull function, 
we will be considering quasi-periodic potentials. 

The mathematical theory of the spectrum of quasi-periodic Schrodinger 
operators is well developed [IPF921 IdlLHlOl . In particular, it is known that 
the spectrum is independent of the space in which it is considered, and, 
more important for us, that the spectrum can be characterized by the ex¬ 
istence of approximate eigenfunctions. In the dynamical interpretation in 
this section, the spectrum corresponds to the Lyapunov exponents of the 
solution nAMB92| . 

In the case of (ITTl) . we can study the Lyapunov spectra for any orbit using 
the geometric constraints (1^ . 

Proposition 4. Let {pn,qn) be an orbit of the mapping given by (ITTl) . As¬ 
sume that it is an orbit in the full measure set that Osledets Theorem applies. 
Then, d—\ Lyapunov exponents are zero. Also, the sum of all the Lyapunov 
exponents is zero. 

Proof. Consider F, the lift of the map F in (1381) . 

Let 5 be a vector perpendicular to a. It is a simple computation to show 
that: 


F{Mq^+f) = + s. 

Then it is clear that the d - \ vectors in the directions perpendicular to s 
do not grow. 

The fact that the sum of the Lyapunov exponents for orbits of a volume 
preserving map is zero is well known since the sums of the Lyapunov ex¬ 
ponents is the rate of growth of the determinant of iterates of the map. □ 

Of course, the dynamical system (1371) is straightforward to implement 
numerically and allows study of statistical properties of depinning. 
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